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ABSTRACT. Changes in ice-sheet surface elevation are caused by a combination of ice-dynamic 
imbalance, ablation, temporal variations in accumulation rate, firn compaction and underlying bedrock 
motion. Thus, deriving the rate of ice-sheet mass change from measured surface elevation change 
requires information on the rate of firn compaction and bedrock motion, which do not involve changes 
in mass, and requires an appropriate firn density to associate with elevation changes induced by recent 
accumulation rate variability. We use a 25 year record of surface temperature and a parameterization 
for accumulation change as a function of temperature to drive a firn compaction model. We apply this 
formulation to ICESat measurements of surface elevation change at three locations on the Greenland ice 
sheet in order to separate the accumulation-driven changes from the ice-dynamic/ablation-driven 
changes, and thus to derive the corresponding mass change. Our calculated densities for the 
accumulation-driven changes range from 410 to 610kgm -3 , which along with 900 kg m -3 for the 
dynamic/ablation-driven changes gives average densities ranging from 680 to 790 kg m -3 . We show that 
using an average (or 'effective') density to convert elevation change to mass change is not valid where 
the accumulation and the dynamic elevation changes are of opposite sign. 


INTRODUCTION 

Knowledge of the mass balance of polar ice sheets is essential 
for understanding sea-level change. The recent assessment 
from the Intergovernmental Panel on Climate Change Fourth 
Assessment Report (IPCC, 2007) summarizes estimates of the 
mass change from both the Greenland and Antarctic ice 
sheets from 1993 to 2003 in the ranges from in balance 
(OGta -1 ) to a loss of 300 Gt a -1 , equivalent to a rate of sea- 
level rise of 0-0.8 mm a -1 . Other results from satellite and 
aircraft measurements of surface elevation change (d/T/df), 
satellite measurements of changes in gravity, and the mass 
input-output method are reviewed (Alley and others, 2007; 
Shepherd and Wingham, 2007) and tabulated in Dahl-Jensen 
and others (2009). A particular issue of concern for deriving 
mass changes (d/VI/df) from observed dH/dt has been the 
appropriate density to use, because dH/dt is in general a 
combination of changes in firn thickness and solid-ice 
thickness, both of which have different densities. Recently 
Zwally and others (201 1) reported a net mass loss from the 
Greenland ice sheet of 1 71 ±4 Gta -1 for the period 2003-07 
and a loss of 7 ± 3 Gta -1 for the period 1 992-2002 based on 
dH/dt observations over these time periods. The methods for 
these calculations are described in more detail in this paper. 

The relation between dH/dt and dA4/df is the result of 
multiple processes occurring throughout the ice-sheet 
column. For example, variations in surface air temperature 
cause a change in the rate of compaction in the upper firn 
layers, which causes a surface elevation change, but this 
temperature-driven change does not involve any mass 
change (Zwally and Li, 2002). Arthern and Wingham 
(1998) evaluated the impact of changes in accumulation 
rate, temperature and surface snow density by using a firn 
compaction model. Their results suggested that accumu- 
lation-rate induced changes were the most significant, and 
that temperature effects were relatively small. In contrast, 
the firn compaction model of Zwally and Li (2002) showed a 


larger temperature sensitivity due to their introduction of a 
temperature-dependent activation energy and rate coeffi- 
cient based on laboratory experiments. Their results showed 
that the seasonal variation in surface temperature caused 
seasonal variations in surface elevation, as well as longer- 
term changes from year-to-year trends in temperature. This 
strong dependence of firn compaction rate on temperature 
has since been supported by field observations of compac- 
tion rates (Arthern and others, 2010). 

Prior to Zwally and others (2011), a constant effective 
density, p e ff, with values between 300 and 900 kg m -3 was 
generally used to convert dH/dt to dA4/df (e.g. 350 in Davis 
and others, 2005; 300 for the 'interior' and 900 to 'seaward' 
in Thomas and others, 2006; 350 and 917 in Wingham and 
others, 2006). FHowever, Figure 1 shows a case in which 
using a constant density is clearly not valid. Consider a 
location where the accumulation rate, A(t), is increasing in 
time while the ice sheet is dynamically thinning (i.e. the firn 
thickness change d/ fim /df>0, and the underlying ice 
thickness change d/ ice /df<0, with example values of firn 
and ice thickness changes and associated densities given in 
Figure 1. The net mass change (per unit area and time) is 
dA4/df=200-900=-700kgrrr 2 a -1 , and the net thickness 
change is d//dt = 0.4-1.0 = -0.6ma _1 . Thus, the effective 
density using this dl/dt would be p e ff=(dA4/df)/(d//dt) = 
-700/-0.60= 1 1 67 kg m~ 3 , which is greater than the density 
of glacier ice, pi = 900kgm -3 . 

Our example clearly shows that the assumption that some 
value of p e ff between 300 and 900 kg m -3 can be chosen is 
not always valid. In fact, it will be invalid wherever there is a 
combination of accumulation-driven thickening and dynam- 
ic thinning, or vice versa. Furthermore, where d/f irn /df 
= d/ ice /df, the corresponding value of p e ff will become 
infinite. This example demonstrates the necessity of esti- 
mating changes in firn thickness and ice thickness sepa- 
rately, and of using the appropriate densities for each. 
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p a = 500 kg nr 3 


p, = 900 kg rrr 3 


Fig. 1 . Example of upper firn and deep ice thickness changes and 
associated mass changes, where 200kgrrr 2 a~' of firn is added to 
the firn-ice column by an accumulation increase and 900 kg rrT 2 a -1 
is subtracted by dynamic thinning. In this case, the conventionally 
defined effective density, p e ff = (dA4/df)/(d//df), is 1167kgmT 3 , 
which would be unrealistically greater than the density of ice. 


In this paper, we present a firn compaction model for 
calculating changes in the rate of firn compaction caused by 
changes in both accumulation rate and temperature. The 
model is used to separate the accumulation-driven and 
ablation-dynamic-driven changes in surface elevation, and 
to calculate the appropriate firn density for accumulation- 
driven changes. It incorporates the processes of surface 
melting, percolation and refreezing. It also includes a critical 
improvement to accommodate a time-dependent accumu- 
lation rate, A{t). We apply the method to three selected field 
locations in Greenland. 


THE FIRN COMPACTION MODEL 

As described in Zwally and Li (2002), the dH/dt of an ice 
sheet is composed of several vertical-velocity components: 
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Here accumulation rate, A(t), is the component that raises 
the surface at a rate of A(f)/p sf , where p s f is the density of the 
surface snow. In general, p s f has a value of ~330kgrrf 5 in 
dry snow regions (Paterson, 1994). In areas with surface 
melt, the meltwater percolates downward in the firn layers 
and the density and thickness of firn layers are modified by 
the meltwater and refreezing (Li and others, 2007; Reeh, 
2008). The melting and refreezing process has been 
incorporated in our model and is described in detail by Li 
and others (2007). A b (f) is the ablation rate occurring only in 
the ablation zone, V\ ce is the vertical velocity of the ice at the 
firn/ice transition, and dB/dt is the vertical motion of the 
underlying bedrock. V fc (t) is the velocity of firn compaction 
at the surface, which is the integral of the compressive 
displacement of the firn layers over the length of the firn 
column. Following the normal usage that surface elevation, 
H, is positive upwards and depth, z, is positive downwards, 
dH/dt, A and dB/dt are positive upwards and V fc , A b and V ice 
are positive downwards in Equation (1). According to 
the mass conservation equation (Li and Zwally, 2002), the 
velocity of firn compaction, V fc , at depth z is given by the 


density, p(z), and compaction rate, dp(z)/dt, as 


V fc (z, t) = 


,z 1 

'dp(z)' 

L, p{z) 

dt 


d z, 
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where Zj is the firn-ice transition depth at which the density 
is ~900 kgrrT 3 . 

Essential to the compaction model is the constitutive 
relation between the compaction rate, dp/dt, and the 
physical variables such as accumulation rate, A(t), and firn 
temperature, T(f), that drive the density change. This relation 
is given by 

^El=K[T(z)]A[ Pi -p(z)\. (3) 

Equation (3) is a semi-empirical relation modified from that 
initially proposed by Herron and Langway (1980) for the 
steady-state case, and A = A(f) is the constant accumulation 
rate in their relation. As discussed by Zwally and Li (2002), 
Equation (3) is based on the idea that the compaction rate at 
depth z is determined by changes in overburden pressure, firn 
temperature and density. However, Equation (3) is not time- 
dependent with respect to A{t) at the surface. Therefore A in 
Equation (3) cannot simply be replaced by A(t) at z=0. For 
example, when A{t) at the surface of the ice sheet is zero, the 
compaction of deeper firn layers continues, but Equation (3) 
would indicate that the compaction rate is zero at all depths. 
Since A at least in part reflects the changes of overburden 
pressure, we replace A by the integral from the surface to 
depth zof A(f-Af)/Afto represent the average change of the 
overburden pressure at depth z, where Af is the time taken for 
the layer to propagate from the surface to depth z. 

The temperature dependence of the compaction rate, 
dp(z)/df, in Equation (3) has been conventionally taken as 
following an Arrhenius relation, i.e. 

K(T) = K 0 e~ E/{RT \ (4) 


where K(T) is the rate factor due to the temperature T in 
kelvin. K 0 is a constant, E is the activation energy and R is a 
gas constant. Experimental results of grain growth and ice 
creep show that the temperature dependence of K{T) is more 
sensitive than Equation (4) under constant K 0 and E, 
indicating that E is actually a function of temperature, E(T). 
Its value increases with temperature (]acka and Li, 1994). 
Since both grain growth and ice creep are involved in firn 
compaction, Zwally and Li (2002) applied the experimental 
results to their firn compaction model. They indicated that 
an increasing value of E{T) will decrease K(T) if K 0 is kept 
constant, which is contrary to the experimental data. 
Therefore K 0 must also be a function of temperature, K 0 {T). 
They modified Equation (4) for firn compaction by introdu- 
cing a temperature-dependent activation energy, f(T), and a 
rate constant, K 0 (T), with an adjustable parameter, (3, i.e. 

K(T) = PK C (T). (5) 

The second term on the right-hand side in Equation (5) is the 
grain growth rate, 

K C (T) = K 0C (T)e~ E W {RT \ (6) 


given by the experiments (Fig. 2a). 

An error in the text of the Zwally and Li (2002) paper 
caused confusion in the application of the theory. By using a 
best fit to the experimental data given by Jacka and Li 
(1994), Zwally and Li (2002) derived the empirical expres- 
sions of grain growth rate, K G (T), and the activation energy, 
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7(7, as functions of temperature, T (Fig. 2a and b). Both 
functions were initially given by Zwally and Li (2002, fig. 3a 
and b). However, K G was mislabeled as K 0 G . Rearranging 
Equation (6) we have 

Kog{T) = K G (T)e E ™ RT l (7) 

Figure 2c gives the empirical form of K 0G {T) as a function of 
temperature by substituting the expressions of E(T) and K C (T) 
shown in Figure 2a and b into Equation (6), together with the 
function of E(T). As shown by Figure 2c, K 0G increases with 
temperature apparently faster than E(T), indicating the 
stronger temperature dependence of K 0C {T) compared to 
E(T). Combining Equations (4-6), we then have 

K(T) = /38. 36(273. 2 - T)“ 2 061 , (8) 

where /3 is an adjustable parameter determined by fitting 
modeled density profiles to field measurements. Equation (8) 
accounts for the temperature dependence of grain growth 
and ice creep rates on firn compaction. It leads to much 
higher temperature sensitivity than Equation (4) under 
constant values of K 0 and E (e.g. Herron and Langway, 
1980), particularly for temperatures higher than -10°C. 
Equations (1-8) are coupled with a one-dimensional heat- 
conducting equation and solved using the multilayer system 
described by Zwally and Li (2002). 

The empirical parameter, (3, in Equation (5) is used to 
calibrate modeled density profiles to field measurements. In 
our previous version of the model, we presented (3 as a 
function of the annual mean temperature based on field 
density profiles from Greenland (Li and others, 2003). 
Helsen and others (2008) extended this relation by using 
41 observed pore close-off depths (depth of the 830 kg rrT 3 
density) from Antarctica where firn temperatures are much 
colder. They found a slightly different relation between f3 
and the annual mean temperature. However, for tempera- 
tures higher than about -30°C (e.g. Greenland), the two 
relations are similar. To further improve the calibration, we 
use two critical points at densities of 550 and 830 kg rrT 3 as 
the control, and tune the value of (3 to force the modeled 
ages at these two densities to match those given by Herron 
and Langway (1980) that were well constrained by the field 
density measurements. Although the previous calibrations 
(Li and others, 2003; Helsen and others, 2008) showed that 
f3 is a function of annual mean temperature only, our tests 
indicate that the accumulation rate also had an influence, 
similar to the Herron and Langway (1980) density-age 
relation for which both temperature and accumulation rate 
were involved. Our present values of (3 as a function of 
annual mean temperature, T m (°C ), and long-term accumu- 
lation rate, (A) (ma -1 ), for Greenland are 

13 = 13 1 

= -9.788 + 8.996(A) -0.61 657 m , (9) 

p < 550 kg m -3 

(3=P2 

= /?1/(— 2.01 78 + 8.4043(A) - 0.0932 T m ), (10) 
p > 550 kg rrT 5 . 

Equations (9) and (10) produce modeled density-depth (age) 
profiles that are in agreement with those from Herron and 
Langway (1980) within an error <±1%. 



Fig. 2. (a) The temperature dependence of grain growth rate, K G ; 
(b) activation energy, E(7), for grain growth; and (c) the derived rate 
constant, K 0G , for grain growth as a function of temperature, 
according to Equation (6) (shown together with E(7) for com- 
parison). Note that the empirical functions in (a) and (b) were 
initially given by Zwally and Li (2002, figs 3b and a). However, K G 
was mislabeled as K 0G in their figure 3b. 


ELEVATION CHANGE COMPONENTS FROM 
ALTIMETRY dH/df 

As described in Equation (1), the surface elevation change, 
d/T/df, is due to a combination of vertical velocity com- 
ponents from different processes. For non-steady state, these 
components can be represented by perturbations from 
steady state, assuming that d6/df is constant during the 
measurements: 

d H d H a dC AT d H d d H b d B 

—T~ = — ; 1 7 — H j— 4 r~ d , 11) 

df df df df df df 

where dH a /df=(A(f)-(A))/p sf is the direct change caused by 
A(f), dC A T/df=-(Vf c (f)-(\/f c )) is the firn compaction-rate 
caused elevation change driven by both A(f) and T(t), 
d/Tb/df=-(A b (f)-(Ab))/pi is driven by changes in the ab- 
lation rate, and dH d /dt=-{V, ce {t)-(V lce )) is driven by 
dynamic changes in the ice flow relative to (A). The () 
symbol indicates long-term averages of the various com- 
ponents, obtained during our model spin-up to a steady 
state. We consider that ablation is zero in the accumulation 
zone where firn exists. In the ablation zone where there is no 
firn, accumulation is zero. In general, dC AT /df depends on 
the history of both A{t) and 7(f) as their effects propagate 
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into the firn. We separate their effects by assuming 

cICat _ dC A dCr 

df df df 

where dC A /df and dC T /df are changes driven by A(t) and 
T(f), respectively. The total elevation change from both A(t) 
and T(t) perturbations becomes 

dH a C AT _ d/T a dC A dCj _ d/-/ a CA dC T 

df — df df df df df 

where d/T^/df is the total change caused by A(t), including 
both the direct change and that associated with the 
compaction rate. Equation (11) then gives 

d l_dH dC T dB dH a CA dH bd 

df - df df df df df ' 1 j 

where d//df is the net thickness change in the firn/ice 
column, defined by treating dB/df and dCT/df as corrections 
to dH/dt. dHbd/df is the combined ablation and dynamic 
term (as dH bd /dt=dH d /dt+ dH^/dt). In the ablation zone 
where there is no firn, dT/ bd /df is the mixed term of both 
dynamics and melting (ablation). In the accumulation zone, 
d/T bd /df is only the dynamic term. The two terms on the 
right-hand side in Equation (12) involve change in mass. To 
obtain d/F^/df, our compaction model first calculates 
d/T a CAT /df using both T(f) and /4(f), and then calculates 
dCj/df using T(f) with constant {A). The dhP CA /dt is then 
given by 

dH a CA _ dH a cAT dCr 

““ d T~ df "d T' 


DERIVING DENSITY, p a , FOR dH a Jdt 

The density, p a , associated with d/T a CA /df shown in Equa- 
tions (13) and (14) is the average density caused by 
perturbations to the accumulation rate over a specified time 
period. Considering the surface elevation change driven by 
variations in accumulation rate /4(f) only, i.e. dH/dt = 
d/T a cA/df, the corresponding mass change, AM a (f), over 
Af=f-f 0 is given by 

AMa(f) = [' [/4(f) - (A)] df. 

Jto 

Here /4(f) and (A) are in units of kgm _2 a _1 . AM a (f) is in 
kgrrT 2 . This mass change causes a surface elevation change 
of 


A Yf 


ca(0 = f 

Jto 


dll' 


CA 


df 


df. 


Thus, 


,,5> 


Equation (15) defines the density, p a , that is obtained using 
the firn compaction model. It represents the density of the 
firn that has been added (removed) as a result of an increase 
(decrease) in A{t) relative to the density profile for the long- 
term average {A}. However, the corresponding change in the 
amount of firn and the compaction rate is distributed 
throughout the firn column, so there is no associated 
physical layer in the firn with a density p a and thickness 
d/T a CA /df. 


i.e. the total accumulation-driven change in surface eleva- 
tion is obtained by subtracting the temperature-driven 
compaction from that caused directly by the change in 
accumulation rate and by the change in compaction rate 
driven by both accumulation and temperature variations. 


MASS CHANGE FROM dH/dt COMPONENTS 

Applying appropriate densities to Equation (12) on the right- 
hand side, the mass change rate, dM/df, of the ice sheet is 
given by 


dM 

df 


(kg m 2 



— — -j— — — T" 900 ~t~ • 
df df 


(13) 


The appropriate density associated with the dynamic and 
ablation term, d/T bd /df (ma -1 ), is the density of glacier ice 
(900 kgrrT 2 ). The density, p a (kgrrT 3 ), associated with the 
shorter-term change of dHWdf (ma -1 ) is discussed in the 
next section. In terms of d//df (ma -1 ), Equation (13) can be 
written as 


dM 

"d T 


900d / 

“dF 


(900 - p a ) 


dH a CA 

df 


(14) 


Equation (1 4) is used to calculate the mass change of the ice 
sheet. It requires the altimetry measurements of dH/dt, 
surface temperature T(t) and accumulation rate A(f) histories 
as inputs to the compaction model to obtain dC T /df, 
d/T a CA /df and the density, p a . Equation (14) indicates that 
using dl/dt with the density of ice to derive dM/df is an 
approximation that neglects the effects of recent changes in 
accumulation rate on the rate of firn compaction. 


RESULTS AND DISCUSSION 

The appropriate method to derive the mass change, dM/df, is 
described by Equation (14). Here we derive the values at 
three selected locations on the Greenland ice sheet (Fig. 3). 
Site A is in the northern region and is associated with lower 
accumulation rate and temperature. Site B is near the 
Summit of the ice sheet, and site C, in the south, is 
associated with much higher accumulation rate and tem- 
perature. We use monthly surface temperatures, T s (f), for 
January 1982-October 2007 as shown in Figure 4a from the 
Advanced Very High Resolution Radiometer (AVHRR) 
(Comiso, 2003). We parameterize changes in the accumu- 
lation rate, A(f), according to the annual mean surface 
temperature at the rate of 5%K -1 . The IPCC report (IPCC, 
2007) summarized the sensitivity of accumulation rate to 
temperature with a range 4. 7-8. 5% K -1 . Data from Green- 
land ice cores (Clausen and others, 1988) give a range 
4-6% KT 1 based on a sensitivity of S ,8 0 to temperature of 
0.69%o5 la OK -1 (Zwally and Giovinetto, 1997). For 1988- 
2005, this parameterization gives an average A(f) increase of 
0.6% a -1 . This is close to the 0.7% a -1 for 1988-2004 
inferred from a regional-climate and surface mass-balance 
model for the percolation and dry-snow zones of Greenland 
(Box and others, 2006). Although there is a good relation 
between changes in accumulation and temperature for 
annual mean values, the correlation could be poor for short- 
term changes (Kapsner and others, 1995). Therefore, in our 
A(f) parameterization, we use the annual mean tempera- 
tures instead of monthly temperatures to avoid introducing 
seasonal variations in A(f) due to strong seasonal variations 
in temperature. A more sophisticated approach that would 
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Fig. 3. A map showing three selected test locations: A (79.4° N, 
31 9.9° E), B (72.9° N, 321. 2° E) and C (61 .9° N, 31 5.9° E) on the 
Greenland ice sheet. The climatic characteristics for the three 
locations are given in Table 1. 


avoid some of these complications in the sensitivity between 
temperature and accumulation rate would be to use A{t) 
from the climate models (e.g. Helsen and others, 2008). 

Starting from the initial steady-state spin-up, we applied 
the time series of 7" s (f) and A(t) to our compaction model. 
Figure 4c shows the modeled time variations in the three 
components from firn compaction at site B, driven by the 
variations in 7~ s (f) and A{t) shown in Figure 4a and b. The 
profiles show significant interannual variations, especially 
since 2000 due to the enhanced warming in that period. 
Although the input temperatures are monthly mean values, 
the annual cycles in the temperature profile still cause 
seasonal variations in C T (f) and /-^cat ( f) profiles, but with 
smaller amplitudes (~2 cm) compared to the results from 
model runs using hourly or daily timescales (Zwally and Li, 
2002). In contrast, the F7 a CA (f) profile does not have seasonal 
variations, because our parameterization of A{t) does not 
have a seasonal variation. We use best fits to the profiles in 
Figure 4c for site B (and the similar profiles for sites A and C), 
for the period October 2003-October 2007, to obtain the 
mass change rates, dM/dt, and associated components for the 
three sites as summarized in Table 1. The associated dH/dt 
are derived from Ice, Cloud and land Elevation Satellite 
(ICESat) measurements for the same period, and dB/df were 




0 5 10 15 20 25 

Year 


Fig. 4. (a) Mean monthly surface temperature at site B in Greenland 
(Fig. 3) showing the general warming during the last two decades, 
(b) The associated increase in accumulation rate as a function of the 
annual mean temperature using a rate of 5%K _1 . (c) The modeled 
surface height changes in three components (as marked) from firn 
compaction driven by temperature and the accumulation rate 
histories shown in (a) and (b). 


Table 1 . The climatic characteristics and the derived rates of the mass change, d/Vf/df, and associated values of the components of elevation 
change as described by Equation (14), together with other parameters for three selected locations (Fig. 3) on the Greenland ice sheet. 
Calculation of dM/df requires determination of the accumulation-driven component, dH^/df, and its associated density, p a , and the ice- 
dynamic- and ablation-driven component, dHbd/df, with density of 900 kg rrT 3 , as well as correction of the observed dH/dt for temperature- 
driven variations in firn compaction, dC T /df, and bedrock motion, dB/df 


Site 

(A) 

Tn 

dH/df 

dcy/dt 

dB/df 

dl/dt 

d /“Cov'd f 

dHbd/df 

Pa 

Pavg 

Pe ff 

d/Wdf 


m a -1 

°C 

m a -1 

m a~' 

m a -1 

m a 1 

m a -1 

m a -1 

kg m -3 

kg m -3 

kg m 3 

kg m~ 2 a -1 

A 

0.10 

-29.3 

0.068 

-0.016 

0.002 

0.082 

0.019 

0.063 

410 

786 

786 

64.5 

B 

0.22 

-29.7 

0.046 

-0.023 

-0.003 

0.072 

0.036 

0.036 

460 

680 

680 

49.0 

C 

0.90 

-14.0 

-0.083 

-0.051 

0.003 

-0.035 

0.104 

-0.139 

610 

776 

1762 

-61.7 
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averaged from three models as in Zwally and others (2005). 
As shown in Table 1 , the impacts from the temperature- and 
accumulation-driven components, dC T /df and d/-/ a cA/df, are 
of similar magnitudes, but the temperature increase causes a 
surface elevation decrease that is offset by the increase in the 
accumulation rate. At site C, the measured surface elevation 
change, dH/dt, is -0.083 m a -1 and d/T a CA /df is 0.1 04 m a -1 , 
indicating that the obtained mass loss (-61 .7 kg m~ 2 a -1 ) is 
due to the ice-dynamic imbalance and/or ablation. Also, we 
note that for site C the values of dM/dt and d//df from 
Table 1 give a conventionally defined effective density of 
1 762 kg m -3 , which is greater than the density of ice as in the 
example of Figure 1 . 

Once the components required for d/VI/df in Equation (14) 
are known, the average physical density can be calculated 
using 


(Pa 

d// a C A 

+ Pi 

dH bd 

) 

d t 

d t 


dH a CA 

+ 

dH bd 


d t 

d t 


Together with Equation (12), Equation (16) describes the 
relationship between the net firn/ice thickness change, d//df, 
and the total thickness change (|d/T a CA /df| + |dA7bd/df|) 
associated with the derived densities. When d/T^/dt and 
dH bb /dt have the same sign, p avg equals p e h (=dA4/df/d//df) 
and the conventional method (d/Wdf = p e ffd//df ) would give 
an arithmetically valid mass change. However, when 
d/T a CA /df and dH bd /dt have opposite signs (e.g. Fig. 1 and 
Table 1 for site C), the conventional method is not valid. 
Equation (16) shows the importance of separating total 
thickness change into firn thickness change and ice 
thickness change. The derived p avg for the three sites are 
summarized in Table 1, together with the corresponding 
values of p e ^ for comparison. The comparison shows that 
p aV g is a more meaningful parameter than p e ff, but neither 
Pavg nor Peff can be used to convert d//df or dH/dt to d/VI/df. 

Positive values of dH bd /dt for sites A and B as shown in 
Table 1 indicate that the ice flow is less than required to 
balance the long-term average (A), which results in dynamic 
thickening in the same manner as negative d/Tbd/df results in 
dynamic thinning. The distribution of accumulation-driven 
changes and dynamic-/ablation-driven changes of the ice 
sheet, and changes in those parameters with time are 
described by Zwally and others (2011). 

CONCLUSION 

Altimetry-derived surface elevation changes, dH/dt, can be 
divided into the components driven by the variations in the 
following: short-term accumulation rate (i.e. d/T^/df), 
surface temperature (i.e. dC T /df), ice-dynamic imbalance 
and ablation (i.e. dH bd /dt) and bedrock motion (i.e. d6/df). 
Partitioning altimetry-derived surface elevation change 
allows assessment of the impacts from recent climate 
variations versus those from long-term ice-dynamic 
changes, leading to more accurate estimations of ice-sheet 
mass balance. Along with 900 kg rrf 3 for the dynamic-/ 
ablation-driven changes, our calculated densities for the 
accumulation-driven changes range from 41 0 to 61 0 kgm” 3 , 
corresponding to average densities, p avg , ranging from 680 to 
790 kgm" 5 . In the locations where thickness changes in the 
firn and ice layers are in opposite directions, the net 
thickness change as derived from the observed dH/dt cannot 


be transformed into mass change simply by using a constant 
average (or effective) density (i.e. a value between 300 and 
900 kg m -3 ). 
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